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Abstract 

The phase-space approach to classical and quantum systems demands for advanced an- 
alytical tools. Such an approach characterizes the evolution of a physical system through a 
set of variables, reducing to the canonically conjugate variables in the classical limit. It often 
happens that phase-space distributions can be written in terms of quadratic forms involving 
the above quoted variables. A significant analytical tool to treat these problems may come 
from the generalized many-variables Hermite polynomials, defined on quadratic forms in TZ n . 
They form an orthonormal system in many dimensions and seem the natural tool to treat the 
harmonic oscillator dynamics in phase-space. In this contribution we discuss the properties 
of these polynomials and present some applications to physical problems. 


1 Introduction 

Classical special functions play a central role in both pure and applied mathematics. Usually 
conceived for the solution of very specific problems, they gave rise to a far-reaching theory, being 
part of, and frequently motivation for, important general theories. 

The trigonometric functions, for instance, originally introduced to deal with specific problems 
of astronomy and navigation, are the basis for the theory of Fourier series and Fourier integral, 
which have applications to many parts of physics. 

Similarly, Bessel functions firstly appeared in mathematical physics in the 1738 Bernoulli’s 
memoir, containing enunciations of theorems on the oscillations of heavy chains. Then, they 
reappeared in mechanical problems as the vibration of a stretched membrane or the symmetrical 
and unsymmetrical propagation of heat in solid cylinders and spheres, as well as in astronomical 
problems, related for instance to the elliptic motion of a planet about the sun. Presently, Bessel 
functions have a very wide field of applications, from abstract number theory and theoretical 
astronomy to concrete problems of physics and engineering [1], 
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Correspondingly, classical orthogonal polynomials (Jacobi, Legendre, Hermite), also introduced 
in connection with astronomical problems, are now of great importance in mathematical physics, 
approximation theory as well as in the theory of mechanical quadrature. In addition, they find 
significant applications in quantum mechanics on the determination of discrete energy spectra and 
the corresponding wave functions in fundamental problems [2] . Orthogonal polynomials are indeed 
essential tool to deal with the problems of the harmonic oscillator and the motion of particles in 
a central field. Furthermore, the classical orthogonal polynomials of a discrete variable are of 

interest in the theory of difference methods [2]. 

Theory of classical special functions is rather well settled [3]. Recursion relations, addition 
theorems, integral representations, generating functions, asymptotic formulae, differential equa- 
tions are collected in an organic body, which is however continuously refined and enriched by new 
investigations and new theoretical approaches [4]. 

Let us recall for instance the method illustrated in ref. [1], which suggests a generalization of 
the Rodriguez formula for the classical orthogonal polynomials, thus allowing to obtain explicit 
integral representations of all the special functions and to derive their basic properties. Similarly, 
the possibility of framing special functions within the context of group theory [4, 5] revealed a 
powerful tool permitting derivation of new results and a rational classification of old results, as 
well as suggesting to introduce new classes of functions, related to the recently discovered algebraic 
systems, such as the supergroups and the quantum groups [6]. 

Furthermore, we note the description of orthogonal polynomials by their recursion relations, 
which once regarded as eigenvalue equations allow to look at orthogonal polynomials from the 
viewpoint of scattering theory [4, 7]. 

Recently, interest in special functions has greatly increased in connection with the possibility 
of generalizing the well known functions of mathematical physics to more than one variable and/or 
more than one index. In this regard, the generalization amounts to introducing functions with 
properties analogous to those of the one-variable counterpart. Generating functions are usually 
the key-note for many-variable generalizations of special functions. 

The multivariable Bessel functions, for instance, originally introduced by Appell [8] in con- 
nection with the problem of the elliptic motion of planets [9], have revealed a wealth of possible 
applications to physical and/or purely mathematical problems, as the scattering of laser radia- 
tion by free or weakly bounded electrons, the emission of e.m. radiation by relativistic electrons 
passing through magnetic undulators [10] as well as problems related to the queuing theory [11]. 
Also, they proved their relevance in multiphoton emission and absorption processes by quantum 
systems, which are of interest for the investigation of squeezed states, the relevant Hamiltonian 
operator containing indeed powers of the annihilation and creation operators [12]. 

Correspondingly, the multivariable generalization of orthogonal polynomials has attracted a 
great amount of interest. In particular, as to the Hermite polynomials, let us recall that in ref. 
[13] a procedure has been developed, which generalizing that proposed by Gould and Hopper [14] 
allows to define multivariable generalized Hermite polynomials, providing a complete orthonormal 
set in Ci(R n ) space of square sommable functions with n variables. 

The present paper concerns with the classical many-variable functions introduced by Hermite 
[15], whose application within the context of the phase space approach to physical problems is 
suggested. Accordingly, in Sec. 2 a general view on the many-variable Hermite polynomials is 
presented. The possibility of exploiting the developed formalism within the context of the phase 
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space picture, which is becoming the context where many physical problems are naturally framed, 
is investigated in Sec. 3. In particular, the basic model of the harmonic oscillator is considered. 
General considerations and extensions are presented in Sec. 4. 


2 Many-variable Hermite polynomials 


One-variable Hermite polynomials Unix) 1 can be defined by means of the relation 





( 1 ) 


where they appear as the coefficients of the series expansion of the exponential of a quadratic form 
defined on the real domain. 

The above expression can also be recast in the more usual form 


e *t-t a /2 



( 2 ) 


from which, exploiting the series development of the exponential function and appropriately rear- 
ranging the summation, the well known expression of the Hermite polynomials in form of a finite 
sum can be easily drawn: 


H n {x) = n\ 


I"/2| 

E 


1=0 


(-) l x n ~ 21 
l\2 l (n — 2l)\ ’ 


(3) 


where [u] denotes the largest integer < v. 

Also, according to the formula for the Taylor expansion, (2) provides the Rodriguez formula: 


Ti n {x) = 


(— )"e z3/2 


<£_ 

dx n 


-* a /2 


(4) 


Taking the derivatives of both sides of (1) or (2) with respect to x and t, it is easy to infer the 
recursion relations 


Tin+i xTi n nH n ^\ , 

'H'n = rihn-l , 


( 5 ) 


linking the polynomial of order n to the contiguous ones. The prime denotes derivative with 
respect to x. 

It is immediate to get from the above relations the differential equation obeyed by Tin’s: 


( d? d \ 

\d^ +X fa) nn = nHn ’ 


which allows to understood Ti n as eigenfunction of the operator 

^ d 2 d 

V= ~^ +X Tx' 


( 6 ) 


(7) 


'The use of the script 7i n to denote the Hermite polynomial is in order to avoid confusion with the more common 
polynomial H n : H n (x) = 2 n/2 H n (\/2x). 
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n being the corresponding eigenvalue. 

The operator V is not self adjoint; the adjoint operator 


T> + = 


dx 2 X dx 


admits the eigenfunction H+, explicitly given as 


K(x) = H n {x)e~ x 'f 2 


( 8 ) 

( 9 ) 


belonging to the same eigenvalue as H n . Accordingly, the functions (H n , H+) form a biorthogonal 
set in the usual sense that they satisfy the orthogonality relation 


f d xH n {x)'H m {x)e x = n\\f2n6 n m ■ (19) 

J — OO 

Generalizing the relation (1) to involve a bilinear form defined in R n , Hermite introduced many- 
variable functions [15]. Adopting a matrix notation, we can write down 

( 11 ) 


where x and h are elements of the vector space R n : 


x = 




( 12 ) 


the superscript T meaning transpose. 

Accordingly, M. is a real n x n matrix: M = (u<j), i, j = 1, • • • ,n, which is required to be 
symmetric: a*j = aji, not degenerate, i.e. det M > 0 and positive definite: an > 0 , i = 1 ,. . . ,n. 
The expression (11) can be rewritten in the alternative form 


x T M h-±h T M h 
€r =~ =— 


£ 


hT 

mi! * 


h™* 

( x ) 

Tflfi. 


(13) 


which is the n-variable analog of (2). 

In passing, it is worth noticing that the above expression suggests considering the more general 
bilinear form 

<t>{x, h) = x T A h - 7:h T Rh (14) 

mi 

with the matrices A, R being symmetric but in general different from each other. In this respect, 
let us recall that the Grassman Hermite functions have been introduced using the above quadratic 
form, with A and B_ being antisymmetric matrices and x, h anticommuting variables. A further 
extension of (14) have been considered in ref. [6], with the intent of obtaining a class of functions, 
related to quantum groups. 
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In full analogy with the one-variable case, taking the derivative of (11) or (13) with respect to 
the components of both vectors x and h, we get the recursion relations satisfied by ‘H mil ..., mn (x): 

I 5Z ®ijXj J (x) 0‘ijTTlj'Hm m ,- 1 m„ (x) , 

V- 1 / i = 1 

d n 

Q Xi ^ m i m„(£) = 2J my-1 mnfe) j (15) 

with t ranging from 1 to n. 

Since corresponding to any quadratic form 0(x) = fMx we can associate the adjoint form 
defined in terms of the inverse matrix, i.e. V’OQ = f M -1 ?, it is possible to associate with the 
s the adjoint polynomials (x). The definition resembles the relation (11), but 

involves the transformed vectors: 


Explicitly, we have indeed 


£ = M x , k = M h . 


(16) 


e -\ (f-*) T M -1 (£-£) _ e -\( r Mr l i y' 




r£/n 




mi! ’ * ’ m n ! mir 


•>77*n 


fe) 


(17) 


If K is the identity matrix, the two polynomials coincide: W mi ,...,m„ = £ mi Furthermore, 

they turn into the product of n one-variable Hermite polynomials, i.e. 


>fm 1> ...,m n ($) = n Wm,N • (18) 

1=1 

Taking the derivative of (17) with respect to x and k, the following set of recursion relations can 
be obtained: 


1 n 

+ (32) %iQq li .__ i q n — AjjQjQgj q._ (x) , 

j=i 

= •,9i-l,...,9n(^) ) (19) 

with A = det^ and Ay the minor relevant to the element ay. 

The relevance of the polynomials G qi „ n is clarified by the orthogonality relation, which can 

be proved in the form 


f Rn d ^ e m n ($)S 9l ,... l9n (x) . (20) 

For the explicit derivation of the above relation the reader is addressed to ref. [8]. 

The orthogonality relation (20) can be conveniently exploited to express a given n-variable 
function p(x) in form of a series involving H mi mn and Gm u ...,m n - Accordingly, let us put 

pfe) = Y te) , ( 21 ) 
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or 


(22) 


p($) — ^ m i mnSmic^mnfe) > 

mi 

with the coefficients A mi m „, being specialized according to (20) into 


i4mi m " (2 tt ) n / 2 mi! ' “ m n ! fw> ^ ’ 

: • • • — -r / dze"5 £T = £ p(x)7i mi) ... imn (£) . 

! 77in! JR" 


\/A 1 


mi r .,m n 


(27r) n / 2 mj! ’ ” m„ 


(23) 


The explicit values of the entries of the matrix M should be suggested by the specific problem 
under study. 

It is needless to say that the theory of many-variable Hermite polynomials is very rich. How- 
ever, the above considerations represent all the machinery we need for testing the possibility of 
using these functions as basis for the phase-space analysis of physical problems. For a more 
detailed discussion the interested reader is addressed to refs. (8, 15]. 


3 Two-variable Hermite polynomials and phase space 
picture of dynamical problems 

Phase space picture is becoming the unifying language for both classical and quantum mechanics. 
Phase space formalism is indeed basic to the Hamiltonian formulation of classical mechanics. In 
this connection, the evolution of a dynamical system is described by a number n of independent 
coordinate variables and on the same number of canonically conjugate momenta. The cartesian 
space of these 2n coordinates is just the phase-space. 

Correspondingly, phase space picture of quantum mechanics is becoming increasingly popular. 
Although, the concept of phase-space is not compatible with quantum mechanics, q and p being 
noncommuting operators, the Wigner phase-space representation allows to overcome this prob- 
lem, since in this representation both the coordinate and momentum variables are c— numbers. 
Accordingly, it is possible to perform phase— space canonical transformations as in the case of 
classical mechanics, which correspond to unitary transformations in the Schrddinger picture of 
quantum mechanics. 

Phase space concept appears therefore as the unifying context, where classical as well as 
quantum mechanics can be naturally framed, thus suggesting the possibility to transfer concepts 
and methods from quantum to classical mechanics and viceversa. 

Furthermore, as discussed in ref. [16], phase— space picture provides the natural language for 
quantum optics as well, offering a geometrical view to coherent and squeezed states as circles 
and ellipses respectively. In this connection, taking advantages from the symmetry of the rele- 
vant Wigner phase space distribution function it is possible to calculate expectation values and 
transition probabilities for the above quoted states [16]. 

Finally, let us recall that phase space is the context where the dynamics of electron beams 
moving through magnetic channels is studied and the Hamiltonian optics can be conveniently 
reformulated. 
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As already remarked, the paper is aimed at investigating the possibility of using the many- 
variable Hermite polynomials, briefly discussed in the previous section, as analytical tool in the 
phase-space approach to dynamical problems. 

In the quantum framework, one-variable Hermite polynomials are intimately related to the 
harmonic oscillator dynamics. It is therefore natural to analyse, as first step for our investigation, 
the harmonic oscillator dynamics in the phase-space representation. 

Let us consider therefore the quadratic Hamiltonian 

H = \v 2 + \KsW , (24) 

the variable s playing the role of time. It is needless to stress the relevance of the above Hamilto- 
nian as basic model for many physical problems as well as approximation in many of the existing 
theories. 

As already stressed, in quantum mechanics the Hamiltonian (24) rules the evolution of a 
harmonic oscillator of unit mass and time - dependent frequency k(s), q and p being the position 
and momentum operators. In classical mechanics, it describes for instance the betatron motion of 
a charged particle through a magnetic quadrupole, as in the ray optics it governs the propagation 
of an optical ray through a nonhomogeneous medium with a quadratic profile of the refractive 
index. 

It interesting to notice that in the configuration space picture the evolution of the system 
described by the Hamiltonian (24) is analysed within the context of conceptually and formally 
different approaches: the Schrodinger equation for the quantum wave function and the Hamilton 
equation of motion for the canonically conjugate variables q and p. Conversely, the Von Neumann 
equation for the Wigner distribution function and the classical Liouville equation for the phase- 
space distribution are of the same form, so long as the Hamiltonian of the quantum system is 
quadratic. Hence, time evolution of the Wigner function can be obtained directly from the solution 
of the equation of motion of the corresponding classical system. Harmonic oscillator provides a 
unique example in which classical and quantum phenomenology overlap to a large extent. 

Let us approach the problem within the context of classical mechanics. Accordingly, the 
Liouville equation for the phase space distribution function p(q, p; s) is immediately written down 
as 

d f d d 1 

^p(g, p; s) = | -p^ + fc(s)^ > p(q, p; s ) , (25) 

with an assigned initial condition: p(q,p,s) = po(q,p). 

As introduction to the forthcoming discussion, let us recall that an invariant quadratic form 

/ = x t, L(s)x (26) 

can be associated to any dynamics described by quadratic Hamiltonians- in canonical coordinates 
and momenta. In passing, it is worth stressing that within a quantum context the quadratic 
form (26) is reported as the Ermakov- Lewis invariant [17], the vector & containing obviously the 
position and momentum operators, whilst in classical mechanics it is reported as the Courant- 
Snyder invariant [18], firstly introduced in the analysis of electron beam motion through magnetic 

channels. In the above expression, the two component vector £_ = is acted by the real 2x2 
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matrix T, which furthermore is required to be symmetric and unimodular: det£ = 1 - Just to 
share the language of accelerator physics, we refer to X. as the Twiss matrix 2 and report it in the 
form 

1 >(a ff)' W 

The entries a,/?, 7, usually named as Twiss parameters, play an important role in designing 
transport channels. 

The quadratic form (26) can be depicted in the phase space as an ellipse, whose size and 
orientation are determined by the Twiss coefficients. The area of the ellipse, which is just the 
value of the invariant /, is usually denoted in accelerator physics as I — ire, e being named as 
the beam emittance. It plays a crucial role in characterizing the quality and the dynamics of the 
e-beam. In a single particle analysis, a,0 , 7 defines the contour of the particle trajectory, as 
in an ensemble analysis they define the second order momentum of the phase space distribution 
function, thus providing information on its extent and maximum localization. Explicitly, 

ei = = (p 2 ) - (p ) 2 , 

e(3 = a 2 q = (c q 2 ) - (q) 2 , 

ea=-a 2 qp = -((qp)-(q)<p)) , (28) 


the averages being understood on the distribution function. Accordingly, the emittance e can be 
given the further meaning: 

e * = a Vvp-A- ( 29 ) 

Let us stress now that the Liouville equation admits as particular solution the distribution function 


p(Q,P\s) = ^exp j-^x T rxJ , (30) 

shaped indeed in form of a Gaussian and therefore with a maximum localization within the ellipse 
of area e. 

The above considerations suggest to use the functions 'H rnn e~ 1Si - =- as basis in the phase- 
plane, the entries of the matrix T being chosen as the second order momentum of the distribution. 
In other words we can use the above quoted functions to approximate a generic function defined in 
the phase space by means of two- variable Gaussians, further modelled by the Hermite polynomials 
7i m , n , which give a different maximum localization. The process is perfectly similar to that used 
for one-variable functions using the Hermite functions. 

Let us consider therefore generic distribution function p(q,p;s). According to the results of 
Sec. 2, we can express p(q,p\s) in form of a series: 

p{Qi p ; s ) = Y, a ™A s )Pm,n(Q,P) » (31) 

m,n 

2 The quadratic form (26) can also been regarded as the transcription in phase space of the quantum invariant 
the vector x being formed by the expectation values of the position and momentum operators and the matrix T 
being linked to the covariance matrix. 
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with a mifl (s) and p^ n (q,p) being explicitly given as 


a m , n (s) = J dqj dpp(q,p\s)Q mtn (q,p) , 
Pm,n(Q>P) = ~~ ■ 


(32) 


The superscript H specifies that the polynomials W m ,n have been used in the series expansion. 

The alternative expansion in terms of the adjoint polynomials Q m ,n(q,P ) can be also used, thus 
leading to 

p(q,p;s) = J2 &m,np£ . n (<7> P ) , (33) 

m,n 

with 


bm,n(s) = J dq J dpp{q,p\s)H m ,n{q,p) , 


Pm,nO?.P) 


1 


2-Kem\n\ 


Sm,n(?,p)e , 


(34) 


and the superscript G specifying that the adjoint polynomials Gm t n have been used. 

The entries of the Twiss matrix appearing in the above expression can be conveniently chosen, 
according to the previous discussion, as the second-order momenta of the distribution function 
Po(q,p) a t the initial time. Then, the emittance e is just obtained according to (29). 

Inserting the expression (31) into the Liouville equation (25), we get the set of equations for 
the coefficients a mfn (s), namely 


(^pp k(7 gq )T7lTlCL rn _ l,n— 1 ^qp [m^TYl l)Um— 2,n kfl(jl l)Um,n-2] 

+ vnxx m ^ i >n +i — kua m . ? (35) 

where the relations 

d d 

dq^ rn,n ~( m “1“ l)Pm+l,n j = ^)Pm,n+l (36) 

and 


PPm y n Pm,n - 1 Vqpijn + l)P m +l,n &qqi n l)Pm,n+ 1 > 

QPm } n ““ Pm—l } n “h cr pp( m d" l)Pm+l,n + &qp{n + l)pm,n+l (37) 

have been used, obtained from the recursive relations obeyed by (see Appendix A). 

Let us consider for instance the particular case where the quadrupole strength does not change 
along the direction of motion: k{s) = fc, and the initial distribution function po{q y p ) has the form 


P° ( ?,P ) = 2TC exp {-^ r £4 • 

(38) 

The coefficients a m| „(0) are then given as 


^m,n(0) = ^m,0^n,0 • 

(39) 
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(40) 


In addition, assuming that the electron beam is matched to the quadrupole, namely 

&qj) = 0 ) <7pp 0 * 

it is easy to prove that the distribution function does not change during the motion: 

p(Q,P\s) = po(q,p) ■ ( 41 ) 

More general situations require the numerical handling of eq. (35), from which the evolution of 
p(q,p\s) can be inferred according to the expansion (31). 


4 Concluding remarks 


The analysis developed in the previous sections has been aimed at testing the possibility of using 
many-variable Hermite polynomials as analytical tool, within the context of phase-space picture 
to dynamical problems. 

The discussion has been limited, for illustrative purposes, to the harmonic oscillator dynamics. 
However, it can be easily realized that Hamiltonians containing higher order terms, accounting 
for non linear forces, can be treated by means of the same formalism as well. 

Also, within a quantum mechanical context, the formalism developed might be used in con- 
nection with the Von Neumann equation, which, as already noticed, rules the evolution of the 
Wigner distribution function W(q,p;s ) according to [19] 


l w ^ P]S)= {-% + k 




Furthermore, let us say that, although in Sec. 3 we have considered 1-dimensional dynamics, 
problems with more than one degree of freedom can be analysed, as, for instance, the motion of 
electron-beams along magnetic channels with transverse coupling. In that case, the radial and 
vertical motions cannot be separated and the dynamics should be analysed in a 4-dimensional 
phase— space, the relevant elements consisting of the conjugate variables x, p x , y , Py- 
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APPENDIX A 

Two-variable Hermite polynomials 


This Appendix is devoted to discuss with some details the properties of the two-variable 
Hermite polynomials, which are used in Sec. 3 as illustrative of the usefulness of such functions 
within the context of the phase-space formalism. 

Firstly, let us write down the explicit form of the recursion relations, reported in (15) for the 
general case of many- variable polynomials: 


W m+ i,„(<7,p) = ( aq + bp)H mi n(q,p ) - amH m -\, n {q,p) - briH m<n -\{q,p ) , 

^m,n+i (?) p) {bq 4" cp)'bi.m,n{.q i p) brriH.m—\,n{.qiP) X (? > P) , (A.l) 

Q 

~Q^im,n{q ) p) (ITTlhLm— 1 ,n (<7, p) "I" bn'H rrltn — l(<7,p) , 


Q^im,n(q,p) = bmH m -\ <n {q,p) 4- cnH m , n -i(q,p) . 

From the above relations after some algebra the following partied differential equation can be 
deduced 


d 2 

C dq 2 


&_ d 2 

a dp 2 dqdp 


+ A 



A(m + , 


(A.2) 


in some sense reminiscent of eq. (6), which can be recovered in correspondence with M = I. 
Consequently, the polynomials H mtTl can be understood as eigenfunctions of the operator 


d 2 d 2 d 2 
T = ~ C !hP ~ a dfi + 2b &idp + A 


( 9 d\ 

\ q dq +P dp) ’ 


(A3) 


with eigenvalue A(m + n). 

Let us consider now the adjoint polynomials Q m ,n(q,p), for which the general relations (19) 
specialize into 


Gm+l,n(q> P) = qGm, n (Q’P ) “ CmGm-l,n(Q> P) + &n£m,n-l (q, p) , 

Gm,n+l(q! p) pGm,n(q , p) 4" bmG m — l,n(.q>P) 1 (q > p) , (A.4) 

'Q^Gm,n{q, P) TTlGm— l,n(<7> p) > Q^Gm,n{q i p) = TlGm,n— l(<7>P) > 

which provide the differential equation 


c dq 2 a dp 2 + 2b dqd P + A (' g dq + p d P ) Grn,n ~ A(m + n ^ m<n ’ 

the same as for H m>n . 

Finally, let us note that the orthogonality relation (20) in the case we are considering specialize 
as 

J d? J dpW min (9,p)^ r>s (<7,p)exp|--x T M£) = -JLm\n\6m,r6n,3 • (A6) 
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